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We add simple tunnelling effects and ray-splitting into the recent trajectory-based semiclassical 
theory of quantum chaotic transport. We use this to derive the weak-localization correction to 
conductance and the shot-noise for a quantum chaotic cavity (billiard) coupled to n leads via tunnel- 
barriers. We derive results for arbitrary tunnelling rates and arbitrary (positive) Ehrenfest time, te. 
For all Ehrenfest times, we show that the shot-noise is enhanced by the tunnelling, while the weak- 
localization is suppressed. In the opaque barrier limit (small tunnelling rates with large lead widths, 
such that the Drude conductance remains finite), the weak-localization goes to zero linearly with 
the tunnelling rate, while the Fano factor of the shot-noise remains finite but becomes independent 
of the Ehrenfest time. The crossover from RMT behaviour (te = 0) to classical behaviour (te = oo) 
goes exponentially with the ratio of the Ehrenfest time to the paired-paths survival time. The 
paired-paths survival time varies between the dwell time (in the transparent barrier limit) and half 
the dwell time (in the opaque barrier limit). Finally our method enables us to see the physical 
origin of the suppression of weak-localization; it is due to the fact that tunnel-barriers "smear" the 
coherent-backscattering peak over reflection and transmission modes. 

PACS numbers: 73.23.-b, 74.40.+k, 05.45.Mt, 05.45.Pq 



I. INTRODUCTION. 

Semiclassical trajectory-based methods have long been 
applied to quantum chaos^, however recent years have 
seen phenomenal progress in this field. The challenge 
of going beyond the Berry diagonal approximation 3 - has 
finally been overcome^. Some works^&i&iLiLii have 
made progress towards a microscopic foundation for the 
well-established conjecture 1 ^ that hyperbolic chaotic sys- 
tems have properties described by random matrix the- 
ory (RMT). Other works 1 ^^^^ 7 .^^^* 21 .^ have 
explored the crossover from the RMT regime (negligi- 
ble Ehrenfest time) to the classical limit (infinite Ehren- 
fest time). Most of the latter (and some of the for- 
mer) works have dealt with mesoscopic effects in the 
quantum transport of electrons through open quan- 
tum dots with chaotic shape. They have analyzed 
the weak-localization correction to conductance and the 
resulting magneto-conductanc e 5 ' 7 ' 12 : 14 : 16 : 18 : 19 , the shot- 
noise (intrinsic quantum noise) 8 ' 13 ' 17 ' 22 , and conductance 
fluctuation a 20 ' 21 . However to-date all such works on 
quantum transport have assumed perfect coupling to the 
leads. The objective of this work is to consider quan- 
tum transport through a chaotic system that is cou- 
pled to the leads via tunnel-barriers (see Fig. [l}. Such 
systems are experimentally realizable as large ballistic 
quantum dots^ 3 -. We study two properties of quantum 
transport, which highlight different aspects of the wave- 
nature of electrons. The first is the weak-localization 
correction to conductance, it is a contribution to con- 
ductance due to interference between electron- waves. It 
is destroyed by a weak-magnetic field, leading to a finite 
magneto-conductance. The second is the shot-noise, an 
intrinsically quantum part of the fluctuations of a non- 
equilibrium electronic current, it is due to the fact that 
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FIG. 1: A chaotic cavity (billiard) with n leads. The mth 
lead has width W m , and is coupled to the cavity through a 
tunnel-barrier with transmission probability, p m . We assume 
that all the tunnel-barriers are very high and narrow, so that 
this transmission probability is independent of the lead mode 
number. 



different parts of the electron-wave can go to different 
places. 

If one assumes the chaotic system is well-described by 
a random matrix from the appropriate random matrix 
ensemble, then it is known how to calculate transport 
through a chaotic system coupled to leads via tunnel- 
barriers 2 ^ 5 .. Under this assumption it has been shown 
that a cavity with two identical leads (both with TV modes 
and with barriers with tunnelling probability p), has a 
Drudc dimensionlcss conductance g Drudc = pTV/2, while 
the weak-localization correction goes like g wl = — p/4. 
Thus the weak localization correction goes to zero when 
we take p — > 0, even when one keeps the Drude conduc- 
tance constant (by taking TV — ► oo such that pN remains 
constant). The formal nature of the RMT derivations 
have made it hard to give a simple physical explanation 
of why this suppression of weak-localization occurs. In- 
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deed the suppression seems counter-intuitive, since we 
know that in a chaotic (ergodic and mixing) system the 
details of the system dynamics have little effect on the 
transport properties. Thus one would expect that trans- 
port through a lead with such tunnel-barriers should be 
equivalent to transport through a lead with N' = pN 
modes without tunnel-barriers. This argument predicts 
the correct Drude conductance, however it would lead 
one to predict incorrectly that the weak-localization cor- 
rection is —1/4, instead of — p/4. The vanishing of weak- 
localization in the limit of strong tunnelling is intriguing 
because it shows that there can be competition between 
two of the most fundamental aspects of quantum me- 
chanics: tunnelling and interference. 

One problem with the random matrix assumption is 
that it tells us nothing about the cross-over to the classi- 
cal limit. The parameter which controls this crossover is 
the ratio of the Ehrenfest time to a dwell time. The 
Ehrcnfcst time (first introduced into chaotic quantum 
transport in Ref. [l2j) is given by the time it takes for a 
wavepacket which is well-localized in phase-space to gets 
stretched to a classical lcngthscale. Since the wavelength, 
Ap = h/pp, is much smaller than the system size, L, the 
stretching is well approximated by the flow of the classical 
dynamics. In a hyperbolic chaotic system this stretching 
is exponential, and happens at a rate given by the Lya- 
punov exponent, A. It is convenient to define the Ehren- 
fest time as twice the time for a minimal wavepacket 
(spread by (XfL) 1 / 2 in position and (XfL/H)^ 1 ^ 2 in mo- 
mentum) to stretch to a classical scale. Typically in open 
chaotic systems there are two classical scales, the system 
size L, and the lead width W; there is an Ehrcnfcst time 
associated with each scal o 26 ' 27 , 

t§} = A -1 1ii[L/Af], (1) 
r° p = X-Hn[{L/X F )(W/L)% (2) 

The first we call the closed-cavity Ehrcnfcst time since it 
is the only Ehrenfest scale in a closed system, the sec- 
ond we call the open-cavity Ehrcnfcst time, because it is 
associated with the opening of the leads. However we 
emphasis that both lengthscales (L and W), and hence 
both Ehrcnfcst times, are relevant to open cavitie o 17 i 18 . 

On times shorter than the Ehrenfest times the dynam- 
ics of the system are not sufficiently mixing to give RAIT 
results for transport quantities. In this article we address 
for the first time how tunnelling affects transport when a 
significant proportion of the current is carried by paths 
shorter than the Ehrenfest times. To do this we must 
include the ray-splitting caused by such tunnel-barriers 
into the semiclassical trajectory-based method. We then 
investigate the effect of tunnelling on weak-localization 
and intrinsic quantum noise (shot noise) in a system with 
finite Ehrcnfcst time. In addition to this we find that the 
method enables us to clearly see the physical mechanism 
by which tunnelling suppresses the weak-localization ef- 
fects. The method presented here was used in Ref. [28| 
to investigate the effect of dephasing voltage probes on 
weak localization. 



This article is organized as follows. Section ITTl contains 
qualitative explanations of the central results in this arti- 
cle (particularly the suppression of weak-localization and 
Ehrcnfcst-independence of the shot-noise in the opaque 
barrier limit). Section IlIII discusses the dynamics of the 
classical paths which are relevant for the semiclassical 
calculations. Section IIVI sets up the formalism for in- 
cluding tunnelling in the semiclassical trajectory-based 
method. Sections fVl and IVT1 then present calculations of 
the weak-localization correction to conductance and the 
Fano factor of the shot-noise, respectively. 

II. QUALITATIVE DISCUSSION OF RESULTS. 

A. Suppression of weak-localization by the 
smearing of the coherent-backscattering peak. 

As noted i n 18 i 19 for a system without tunnel-barriers, 
the coherent-backscattering approximately doubles the 
weight of all paths returning close and anti-parallel to 
the same path at injection. The approximate doubling 
of the weight applies to all paths that return to a strip 
defined by (y, 8) across the lead, where 

9 - ^ -p F 1 mX(y - y ) cos 9 (3) 

for a path that was initially injected at (j/q, Oq)- This strip 
sits on the stable axis of the classical dynamics, with a 
width in the unstable direction of order hljt-pW)^ 1 . The 
probability for the particle to go anywhere else is slightly 
reduced (by the weak-localization correction which is 
negative and (y, 0)-independent), this ensures that the 
probability for the particle to go somewhere remains one 
(preserving unitarity and conserving current). Since the 
enhancement due to coherent-backscattering only occurs 
for reflection, we see that reflection is slightly enhanced 
and transmission (and hence conductance) is slightly re- 
duced; this is the weak-localization correction to conduc- 
tance. 

However once we introduce tunnel-barriers the situa- 
tion becomes more complicated. There is still the en- 
hancement of paths that return to close and anti-parallel 
to themselves at injection, however there is no longer 
a guarantee that this enhancement contributes to re- 
flection. The enhanced paths can return to the injec- 
tion lead but then be reflected off the tunnel-barrier, 
remaining in the cavity (see Fig. |4j;). These enhanced 
paths will then bounce in the cavity until they eventu- 
ally transmit through a barrier, either contributing to 
transmission (as shown in Fig.H^) or reflection (as shown 
in Fig. |4j: if we set m = mo). We use the term suc- 
cessful coherent-backscattering to refer to the the usual 
coherent-backscattering contribution, shown in Fig. 0)3, 
while using the term failed coherent-backscattering for 
those contributions which behave in the manner shown 
in Fig. gb,d. 

Thus the failed coherent-backscattering gives a positive 
contribution to both reflection and transmission which 
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will tend to off-set the negative one coming from the usual 
weak-localization correction. This contribution will go 
like (1 — p) times the cohcrent-backscattering contribu- 
tion in the absence of tunnel-barriers. This contribution 
will be "smeared" over all leads with a weight for each 
lead equal to the probability of escaping through that 
lead. So if the cavity has two identical leads, half of it 
will contribute to transmission and the other half to re- 
flection. The successful cohercnt-backscattering goes like 
p and only contributes to reflection. 

As we take p — ► we see that the coherent back- 
scattering becomes completely "failed" (the successful 
part goes to zero). So the whole backscattering peak 
is smeared over all leads in exactly the same way as the 
conventional weak-localization correction. Both have a 
weight for any given lead equal to the probability for 
the particle to escape through that lead. However the 
two have opposite signs (which they must to preserve 
unitarity) and thereby cancel. Hence the total weak- 
localization correction (sum of the conventional weak- 
localization correction, Fig. 2^, and the failed coherent- 
backscattcring, Fig. 0t,d) vanishes when we take p — > 
(even if we keep pN finite, so the Drude conductance 
remains finite). 



B. RMT to classical crossover given by ratio of 
Ehrenfest time to paired-paths survival time. 

In all calculations in this article we find that the 
crossover from RMT behaviour (zero Ehrenfest time) 
to classical behaviour (infinite Ehrenfest time) is gov- 
erned by the ratio of the Ehrenfest time to the paired- 
paths survival time, not the single-path survival time (the 
dwell time). The paired-paths survival time, defined in 
Eq. (|13[) , is the time during which a pair of initially iden- 
tical classical paths both remain inside the cavity. Under 
the classical dynamics of a closed cavity a pair of ini- 
tially identical paths remain paired (identical) forever. 
However the presence of the leads means that all paths 
remain a finite time in the cavity, wc call this the dwell 
time or single-path survival time, tdi- Further, if there 
are tunnel-barriers on the leads then the paired-paths 
survival time, td2, is given by the time for one or both of 
two (initially identical) paths to escape. We show that 
this time varies between tdi for transparent barriers to 
t^tdi for opaque barriers. 

The reason why td2 is the relevant timescale for the 
suppression of the crossover is as follows. All paths which 
contribute to the RMT limit must involve an intersection 
at an angle of order (X-p/L) 1 ^ 2 , and the paths must the 
survive in the cavity until they spread to a distance apart 
of order the width of the lead on both sides of the inter- 
section. Within |Tg P on either side of the intersection the 
paths are correlated, and so their joint survival probabil- 
ity decays at the rate TjjJ ■ The shot noise contributions 
divide into those in which the paired-paths survive for a 
time less than Tg P and those which survive a time greater 



than Tg P . Thus we see that these contributions go like 
(1 — exp[— Tg P /rD2]) and exp[— t^ p /ty>2] respectively. 

All weak-localization contributions require that the 
path segments diverge to a distance of order the system 
size on one side of the intersection (so that a closed loop 
can form) and to a distance of order the lead width on the 
other side (when the "legs" can escape). Thus the path 
segments must survive until a time Tg on the loop-side of 
the intersection and Tg P on the leg-side. However within 
Tg P /2 on both sides of the intersection, the paths have 
the paired-paths survival probability (because they are 
closer than the lead width), while elsewhere each path 
segment individually has the single-path escape proba- 
bility. Thus we see that the exponential suppression of 
weak-localization for finite Ehrenfest time goes like the 
survival probability for such paths of length r^ p + Tg 
which is 

cx P [-t° p /t D2 -(t£-t° p )/t d1 ]. (4) 

The first term in the exponent is the Ehrenfest time de- 
pendence, the second term is a classical correction, since 
(Tg — Tg P ) is independent of the particle wavelength. In 
the transparent barrier limit (where td2 = Tdi) the sup- 
pression goes like exp[— Te/tdi] as found in Ref. [18l ]. 
while in the opaque barrier limit we find td2 = Tdi/2 so 
the suppression goes like exp[— (r^ 1 + t-£ P )/tdi]. 

In general the survival probability of such intersecting 
paths of total length t > 2r^ p is 

cxp[-(t-2r° p )/r D1 -r° p /r D2 ], (5) 

since the paths arc paircd-during a time Tg P and unpaired 
at all other times (when unpaired each segment decays 
at a rate of t^ 1 ). In the transparent barrier limit this is 
exp[— (t — Tg P )/rDi] while in the opaque barrier limit it 
is simply exp[— £/tdi]. 

C. Ehrenfest time independence of shot noise in 
the opaque barrier limit. 

In the limit of opaque barriers (small tunnelling rate) 
we find that the shot noise is independent of the Ehren- 
fest time up to first order in p, see Eq. ([75)1. Thus for 
any small p, the magnitude of the difference between the 
noise in the RMT limit (t^P = 0) and the classical limit 
(t^ p = oo ) goes like p 2 , and thus may be beyond experi- 
mental resolution. This is despite the fact that the con- 
tributions in the RMT limit are very different from those 
in the classical limit. Closer inspection shows that those 
contributions which survive in the opaque limit (Fig. [6J3- 
i and Fig. [7j3-g) do have one thing in common; none of 
them have correlations between the paths when enter- 
ing or exiting the cavity. Thus it appears that it is ir- 
relevant whether correlated paths have enough time to 
become uncorrelated inside the cavity (as they do when 
Tg P <C td2, but not when r^ p 3> 102) because the tunnel- 
barriers destroy all correlations between paths (the prob- 
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ability that both paths exit when they hit the tunnel- 
barrier is zero in this limit). It is intriguing that the 
shot noise is so insensitive to the manner in which corre- 
lations between paths are destroyed (either via classical 
dynamics or via tunnelling). 

We note that for a symmetric two-lead system in the 
opaque limit (Ni = N2 and pi = p% -C 1), the shot noise 
is half that for Poissionian noise (Fano factor, F = 1/2). 
Thus even through the transmission probability is small, 
so transmission is rare, the transmission is sub-Poissonian 
process. The Fermionic nature of the electrons induces 
correlations even in this limit. The result that F = 1/2, 
is the same as one would get for two leads coupled by a 
single tunnel-barrier with transmission probability equal 
to 1/2 (with no chaotic system present at all). 



III. CLASSICAL DYNAMICS IN A SYSTEM 
WITH TUNNEL-BARRIERS ON THE LEADS. 

In this article we consider a quantum system whose 
classical limit (wavelength Af —> 0) has a non-zero tun- 
nelling probability at tunnel-barriers (we assume the 
width of the barriers scales with Ap). Thus it is natural to 
include tunnelling at the level of the classical dynamics. 
Here we introduce tunnelling into the classical dynamics 
phenomenologically, we will then see in section IIVI that 
it is the classical limit of the quantum problem we wish 
to solve. We assume that each time a classical path hits 
a barrier there is a probability of p that the path goes 
through the barrier (as if the barrier is absent) and a 
probability of (1 — p) that the particle is specularly re- 
flected (as if the barrier is impenetrable). 

The probablistic nature of this tunnelling makes the 
classical dynamics stochastic. The classical paths, which 
are the solutions of these stochastic dynamics, have prop- 
erties not present in deterministic dynamics, such as bi- 
furcations (ray-splitting) at the tunnel barriers (one part 
transmitting and the other reflecting). To simplify the 
considerations in this article we consider tunnel-barriers 
that have a tunnelling probability which is independent 
of the angle, 9, at which a classical path hits the bar- 
rier. This assumption is justified for tunnel-barriers in 
the limit of large barrier height and small barrier thick- 
ness (see Appendix \E§. This is equivalent to saying we 
consider that the tunnel-barriers have the same tunnel 
probability for all lead modes. 

Let us now consider the dynamics of a classical particle 
in such a system. We define P(Y,Yo;i) such that the 
classical probability for a particle to go from an initial 
position and momentum angle of Yo = (yo,0o) in the 
phase-space of lead mo (a point on the cross-section of 
the lead just to the lead side of the tunnel-barrier) to 
within (5y, 89) of Y = (y, 9) on the cross-section of lead 
m (again just to the lead side of the tunnel-barrier) in a 
time within 5t of t is 

P(Y,Y ;t)6y696t. (6) 



In a system without tunnel-barriers the classical dynam- 
ics are deterministic, and P(Y,Yo;i) has a Dirac 5- 
function with unit weight on classical paths. However 
in a system with tunnel-barriers the classical dynamics 
are stochastic, with each barrier acting as a ray-splitter. 
A classical path which hits a barrier is split into two 
(one which passes through the barrier and one which 
reflects). Thus, P(Y,Yo;i) has a 5- function on each 
classical path which exists (counting all possible ray- 
splittings), however the weight of the Dirac S- function 
is the product of the tunnelling/reflection probabilities 
that that path has acquired each time it has hit a tunnel- 
barrier. Thus the integral of P(Y, Yo;i) over all posi- 
tions/momenta at escape, Y, gives a sum over all paths 
starting from Yo in which each term is weighted by these 
tunnelling/reflection probabilities, hence 

J^dt J dYP(Y,Y ;t) [•••] Yo 

= E PmP mo n m '[l-Pm'p im,) [---h (7) 
7£{Y } 

where dY = dyd9, with y integrated over the cross- 
section on each lead and 9 integrated from — ir/2 to 
7r/2. The sum is over all possible path starting at 
Yo = (yo,Py a )- Lead m is defined as the lead the path 
finally escapes into, and n 7 (m') is the number of times 
that the path 7 reflects off the tunnel-barrier on lead m' 
before escaping. We assume that any quantities in [■•■] 7 
are independent of p m i , in-other- words they are same as 
for the path which would exist if the tunnel-barriers were 
impenetrable for each reflection and absent for each tun- 
nelling of path 7. 

In semiclassics one is typically summing over all paths, 
7, from yo to y with energy, E, rather than all paths 
starting at Yo. Using Eq. f7]) we see that this sum can 
be written in terms of P(Y, Yo;i) as 

E PmPmoIlm' t 1 ~ Pm']"^" } [• ' ' ]<y 

= J^dt £^MoM(Jj£)P(Y,Y -,t) [•■•K (8) 

where the factor of (dy/d9o) comes from the change of 
integration variable from y to 9q. 



A. Ensemble or energy average of classical 
probabilities. 

If we average, (■••), over energy or an ensemble of 
similar chaotic systems, the average of P(Y,Y ;t) will 
be a smooth function. If the system is mixing (when the 
leads are absent) and we perform sufficient averaging, 
we can assume the probability to go to any place in the 
cavity phase space is uniform. From this it immediately 
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follows that the probability to go from a point in lead too 
to anywhere in lead m (for to ^ m ) is 

[°°dt / dY(P(Y,Y ;t)} = Pma x - PmWm 

PmgPm^m /q\ 

L-im 1 Pm' ^m' 

where Y is integrated over the cross-section of lead to. 
To see this we simply note that the probability to en- 
ter the cavity from lead too is p mo , and the probability 
to escape the cavity is one. Thus the probability to es- 
cape into lead to is the ratio of p m W m to the sum of 
Pm'W m i over all leads. Note that in doing this we ignore 
all unsuccessful attempts to escape, thus the particle may 
have hit tunnel-barriers on the leads many times, but 
each time been reflected. The result in Eq. (J9j) can also 
be derived by noting that the probability of hitting the 
tunnel-barrier on lead rrii is W m j '(tqL) where tq is the 
time of flight across the cavity, and then explicitly sum- 
ming all reflections off the barriers. However we feel the 
above logic of ignoring all unsuccessful escape attempts 
gives the results in a direct and simpler manner, thus we 
use it throughout this article. 

Following the same logic that leads to Eq. ((9]), we see 
that the survival probability of a single classical path 
which is inside the cavity is given by the master equation, 

A(t) = -r D1 1 p 1 (t) (io) 

where we define the single path dwell time as 

n 

T^l = (roi)" 1 PmWm- (11) 

m— 1 

Thus the probability to tunnel into the cavity from lead 
mo then escape into lead to at position Y = {y, 9) on 
the phase-space cross-section just to the lead side of the 
tunnel-barrier is 

<P(Y,Y ;i)) = Pm <> Pm ™ s6 ex P H/r D1 ].(12) 

Now we turn to the evolution of a pair of almost iden- 
tical paths. A quick glance at the figures in this article 
show that such pairs of paths (vertically cross-hatched 
regions in Figs.[3][n]) are crucial to the quantities we cal- 
culate. This situation is different from the evolution of 
two very different paths, because here the two paths ex- 
plore the same regions of the cavity's phase-space. Thus 
if one path hits a tunnel-barrier then so will the other 
one. This leads to the definition of "almost identical" ; 
the paths must be a perpendicular distance apart that is 
less than the lead width for a significant period of time 
(at least a few bounces), i.e. their difference in momenta 
is less than pp x W/ L. (In principle two paths with very 
different momenta, will be "almost identical" very close 
to the points where they cross, however we ignore this 
as it will not significantly affect our calculations.) Then 



when the pair of paths hit a tunnel-barrier the probabil- 
ity to survive (probability for both paths to remain in 
the cavity) is (1 — p m ) 2 i this means the probability to 
not survive is p m (2 — p m )- Thus the survival probability 
of paired-paths is given by the master equation 

w) = -r^w) (is) 

where we define the paired-paths survival time as 

n 

t£ = (toL)- 1 p m {2 - p m )W m . (14) 

m— 1 

Note that this is the probability that both paths sur- 
vive. Since < p m < 1 we see that 1 < tdi/td2 < 2, 
with td2 = tdi when all tunnel-barriers are transparent 
(p m = 1 for all m), and rrj2 = i~di/2 in the limit of 
opaque barriers (p m — > for all to). Unfortunately we 
cannot write down an equation for the evolution of pairs 
of paths which is as simple as Eq. (|12p . because even 
under averaging the position/momentum of the second 
path in the pair are deterministically given by its initial 
position/momentum relative to the first path. Assuming 
the system has locally uniform hyperbolic dynamics (on 
lengthscales up to W <C L), then the dynamics of the 
second path in coordinates perpendicular to the first are 

£ ( "'I 1 )= M ( r± l L ) ™ 

dt \ P±/PF J \P±/PF J 

where M. is a two-by-two matrix. For Hamiltonian dy- 
namics M. has eigenvalues ±A where A is the Lyapunov 
exponent, and so the dynamics are area preserving. We 
will assume for simplicity that the eigenvector of A4 with 
eigenvalue ±A is along the axis defined by = ±m\r±. 
We will use this to calculate the relative position of the 
second path in a pair when it is needed in the calculations 
below. 



IV. TRAJECTORY-BASED SEMICLASSICS 
WITH TUNNEL-BARRIERS. 

The semiclassical derivation of the Drude conductance 
has become standardise (here we will broadly follow 
Ref. [I!])- However we wish to introduce tunnelling into 
this derivation, this presents a difficulty since tunnelling 
cannot be described in conventional scmiclassics. To deal 
with this we follow a well-established procedure for deal- 
ing with a semiclassical system (wavelength much less 
than other lengthscales) in which there are isolated re- 
gions where semiclassics fails, see for example Ref. [30| . 
We treat the regions where semiclassics fails in an exact 
manner (or using an appropriate approximation scheme) 
and then couple the propagators in those regions with 
the semiclassical ones in the regions where semiclassics 
works well. 
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A. The energy Green's function and scattering 
matrix elements. 

Before addressing the construction of an energy 
Green's function which includes the tunnelling, let us 
consider the case of a ray inside the cavity hitting the 
tunnel-barrier on one of the leads. A ray being a plane- 
wave multiplied by an envelope-function in the direction 
perpendicular to its motion which is much wider than Af, 
but much narrower than the classical scales, W, L. In this 
case we can treat the scattering of the ray in the same 
manner as a plane-wave hitting the tunnel-barrier on a 
lead of infinite width. The equations for motion parallel 
and perpendicular to the barrier can be solved separately. 
The evolution perpendicular to the barrier is given by the 
solution of a textbook one-dimensional tunnelling prob- 
lem (see Appendix [A|) . while the evolution parallel to the 
barrier is unchanged by the presence of the barrier. The 
solution of the one-dimensional tunnelling problem tells 
us that for any given ray arriving at the tunnel-barrier, 
there will be two rays leaving the tunnel-barrier. One ray 
is that transmitted through the barrier, it has complex 
amplitude t and has the same momentum as the incom- 
ing ray. The other ray is that reflected off the barrier, it 
has complex amplitude r and has its momentum perpen- 
dicular to the barrier reversed (specular reflection) . Thus 
each ray will be split into two each time it encounters a 
tunnel-barrier. Of course the amplitude of the two new 
rays are such that |r| 2 + |i| 2 = l. To fit with the classical 
model in Section Hill we define p = \t\ 2 . 

For a very narrow barrier (which is the only case we 
consider in this article) these two rays are identical to 
the rays one would have if the barrier was both absent 
(transmitted ray) and impenetrable (reflected ray), ex- 
cept that here the amplitude of the rays are multiplied 
by the complex amplitudes t and r respectively. For such 
narrow tunnel-barriers, the amplitudes t and r are inde- 
pendent of the angle, 9, of the incoming ray. 

Now we must do the same for the energy Green's func- 
tion for propagating from ro in the L lead to r in the R 
lead. We do this in a manner similar to Ref. [3(J. We first 
note that at large distances, |r — ro 3> Af, a scmiclas- 
sical Green's function, G(r,r ; E), is well approximated 
by a plane-wave in the vacinity of any given classical 
path. Further for r 7^ ro the differential equations for the 
Green's function and a wavefunction are the same (the 
Schrodinger equation). Thus we treat the semiclassical 
Green's function in the vacinity of a given path which 
touches the tunnel barrier (path a), as a plane- wave or 
ray. We use this as the ingoing boundary condition on 
cross-section A (see Fig. [2]) for the tunnelling/reflection 
problem. The tunnclling-reflection problem is solved 
using the standard method for wavefunctions (see Ap- 
pendix [XJ . The transmitted part gains a complex pref- 
actor, t, and then couples to the lead modes (as in the 
case without tunnel barriers, we assume the leads are 
wide enough that they accept all momentum states). 
The reflected part gains a complex amplitude, r, and 




FIG. 2: Schematic of the calculation of Green's functions in 
the vacinity of a tunnel barrier. The evolution from r' — > r" 
and r' — > r'" is treated exactly, while semiclassics is used 
for the evolution from ro — » r' and r" — > r. We calculate the 
contribution for the part of the incoming Green's function due 
to paths in the vacinity of the classical path a, by matching 
the semiclasical and exact solutions of the Green's function 
equation on the cross-section A, which is just to the left of 
the tunnel-barrier. Just to the right of the tunnel-barrier, on 
cross-section B, we couple to the lead modes. 



is a plane-wave with its momentum perpendicular to the 
barrier reversed compared with the incoming plane- wave. 
Treating this outgoing plane wave, at cross-section A, as 
a boundary-condition on the semiclassical evolution, we 
see it couples to exactly the same paths as if the barrier 
was impenetrable. Hence every tunnel-barrier couples 
each incoming path to two outgoing paths (one as if the 
barrier was absent and one as if the barrier were impene- 
trable). If the weight on the incoming path is B a then the 
weight on the two outgoing paths will be tB a and rB a 
for transmisson and reflection at the barrier, respectively. 

Hence when we write the semiclassical Green's func- 
tion inside the cavity in the usual way, it is a sum over 
classical paths inside the cavity which can either transmit 
or reflect at each barrier. The properties of the classical 
path (action, S-y, Maslov index, /x 7 , and classical stabil- 
ity) are the same as for the classical path that would 
exist if the barrier were absent for each transmission and 
impenetrable for each reflection. The full energy Green's 
function, including the tunnel-barriers, then takes the 
form ^ 7 B 1 exp[iS~ ( /h + iyU 7 /2], where B 1 is the square- 
root of the stability of the path multiplied by a complex 
factor of t m i and r m i for each transmission and reflection 
at barrier m'. Using Ref. [2!| to go from this Green's 
function to the formula for S mmo; ji, the scattering ma- 
trix element to go from mode i on lead mo to mode j on 
lead m, we get 

§mm ;j. = — (27rift) ~ 1 /2 / dy / d^y^ A y 
Jl Jr ~ 

x(i|2/}(j/oH)exp[i5 7 /?i + i7r/i 7 /2] , (16) 

where \i) is the transverse wavefunction of the ith lead 
mode on lead m! . The sum is over all paths 7 (with clas- 
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sical action Sj and Maslov index /i 7 ) which start at yo on 
the cross-section of the injection (mo) lead, tunnel into 
the cavity, bounce many times (including reflecting off 
the tunnel-barriers on the leads) before tunnelling into 
lead m at point y on its cross-section. The complex am- 
plitude A-y is 



A, 



dpy, 

dy 



1/2 



i n 7 (m 



(17) 



where path 7 starts from Y = (yo,Py ) 011 ^ead uiq and 
tunnels into the cavity, it then bounces inside the cavity, 
finally it transmits through barrier m to end at Y = 
(y,p y ) in lead m. We define n 7 (m') as the number of 
times that the path 7 reflects off the tunnel-barrier on 
lead m! before escaping. The differential {dp Vo /dy) 1 is 
the rate of change of initial momentum, p yo , with final 
position, y, for an unchanged set of transmissions and 
reflections at the barriers. 



B. The Landauer-Biittiker formula and the Drude 
conductance. 

Inserting Eq. (|16p into the Landauer-Biittiker formula 
for the conductance 



9mniQ 



— Tr [S TOm §mm ] j 



(18) 



one gets a double sum over paths, 7I and 72 and over lead 
modes, \n) and |m). We make a scmiclassical approxima- 
tion (see Appendix [B]) that ^2 n (y'\n)(n\y) ~ 8{y' — y). 
The conductance is then given by a double sum over 
paths which both go from yo on lead mo to y on lead 



QmrriQ 



(2nh)- 1 [dy [dyY] A 7l A* l2 e iSS / h . 
l Jn 7 i )7 2 



(19) 



where the action difference, 5S = iS 7 i — S 7 2. 

We averaged the conductance over the energy or the 
cavity shape. For most {71,72} the phase of a given 
contribution, SS/h, will oscillate wildly with these vari- 
ations, so the contribution averages to zero. The only 
contributions that will survive are those in which 7I and 
72 are correlated in such a manner that 6S remains fixed 
when we vary the energy or the cavity shape. The most 
obvious such contributions^ are diagonal ones (7I = 72), 
and we now show that they give the Drude conductance. 

For these diagonal contributions, we note that 
(dp ya /dy) 1 = pf cos 9o{d6o/dy) 1 and then use Eq. ([5]) 
with p m i = \t m '\ = 1 — l^m/l 2 , to write the sum over all 
paths 7 from yo to y as 



,7T/2 

dt / d9 d6 

J-ir/2 



x?(Y,Y ;t)[-]Y . (20) 



where for compactness we define P(Y,Yo;i) = 
Pf cos 8 Q x P(Y,Y ;t), with ppcos(?o being the initial 
momentum along the injection lead. Using Eq. ^ we 
see that 

/E^v V *U PmoPmPF COS O COS 9 

P (Y,Y *)) = — ^ =— exp -t/T D i .(21) 

2 (Lm' Pm'W m ')TDl 

Using Eqs. ([20]) and ([21]) one find the Drude conductance, 

D (-i _ \j\t r 1 Pmo Nm a Pm N m . . 

The first term above is for reflection (<5 m?no is a Kronecker 
delta-function), it represents contributions which are re- 
flected off the barrier on the injection lead without ever 
entering the cavity. 



V. WEAK-LOCALIZATION CORRECTION TO 
CONDUCTANCE FOR A CHAOTIC SYSTEM 
WITH TUNNEL-BARRIERS. 

We identify four contributions to the weak-localization 
correction to the conductance of a system coupled to 
leads via tunnel-barriers, they are shown in Fig. 3J 



A. Usual weak-localization. 

Here we consider the usual contribution to weak lo- 
calization (wlO), the first contribution shown in Fig. 2] 
We call it the usual contribution because it is the only 
weak-localization contribution to transmission (m ^ mo) 
when the tunnel-barriers are absent (p m = 1 for all m). 
In the usual contribution, the paths are paired almost 
everywhere except in the vicinity of an encounter—. At 
an encounter, shown in detail in Fig. [3] one of the paths 
intersects itself, while the other one avoids the crossing. 
Ref. [1] showed that every self-intersecting path with a 
small crossing angle e has a partner which avoids the 
crossing. As a result the two paths travel along the 
loop that they form in opposite directions. However the 
two paths are always close enough to each other that 
they have the same stability, hence 72 ^71^72 ~~ * 
^ 7l |^4 7 i| 2 . To evaluate the usual weak-localization 
contribution (wlO), we perform a calculation similar to 
Ref. [3; this broadly follows Refs. [Hi] while adding 
the crucial fact that close to the encounter (the cross- 
hatched region in Fig. 2} the paths have a different sur- 
vival probability from elsewhere. For a system with 
the tunnel-barriers the survival probability is given by 
Eq. (|13p inside the cross-hatched region, and given by 
Eq. ([TO]) elsewhere. 

As in Ref. we only consider those contributions 
where the paths are uncorrelated when they hit the lead 
(by which we mean a collision with a lead where at least 
one pair of paths tunnels out of the system). Those 
contributions in which the paths are paired when one 
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encounter encounter 



FIG. 3: In all figures of contributions in this article we mark the regions where paths are correlated with vertical cross-hatching. 
The encounter is at the centre of those regions. The diagonal approximation is good everywhere except at the encounter, thus 
for simplicity we show only two of the four paths (or path segments) in all other figures in this article. Here we show what the 
cross-hatched regions in the other figures mean in terms of all four paths that enter them. In the cross-hatched region the two 
paths marked by solid lines are paired with each other. This region is macroscopic in size, the length of the region in time is t^ p 
(which we assume to be of similar magnitude to the classical timescales such as tdi,td2), and the distance between the paths 
at either end of this region is of order the lead width. In contrast the encounter (where both paths marked by dashed lines 
swap from one path marked by a solid line to the other) is microscopic in size, the distance between the non-crossing paths is 
typically (Af£) 1//2 which goes to zero in the classical limit, the time which the encounter takes is a few times the inverse of the 
Lyapunov exponent. 



of them tunnels out of the system are included in the 
coherent-backscattering contributions. Uncorrelated es- 
cape requires a minimal time 7V(e)/2 between encounter 
and escape, where^ 

T w (e)^X- 1 ln[e- 2 (W/L)% (23) 

for a small crossing angle e (see Fig. [3]) , in a system with 
Lypunov exponent A, system size L, and lead width W. 
This is because this is the time for the perpendicular 
distance between the paths to become larger than the 
width of the leads. Only then will the two paths escape in 
an uncorrelated manner, typically at completely different 
times, with completely different momenta (and possibly 
through different leads). 

To calculate the contribution of the sum over all paths 
of the form sketched in Fig. [4j we note that the action 
difference, 5S, is~& 

5S wi = E F e 2 /X. (24) 

This assumes the dynamics in the vacinity of the en- 
counter are time-reversal symmetric, this is the case for 
any applied magnetic field which is weak enough not to 
affect the classical dynamics. We write the probability 
to go from Yo to Y in time i, in terms of the product 
of two probabilities. The first is probability to go from 
Yo to a point on the energy surface inside the cavity 
Ri = (ri , 0i ) (where 0i defines the direction of the mo- 
mentum) in time t\. The second is the probability to go 
from Ri to Y in time t — t±. When one integrates this 
product over Ri on the energy surface C, one gets the 
probability to go from Yo to Y in time t. Thus we can 
write the quantity P, introduced above, as 

P(Y,Y ;t) = / dR 2 dRiP(Y,R 2 ;t-i 2 ) 
Jc 

x P(R 2 ,R 1 ;t 2 -i 1 )P(R 1 ,Y ;ti), (25) 

where P(R 2 , Ri; t) is the probability density to go from 
Ri to R 2 in time t, but P(Ri,Y ;i) is a probability 
density multiplied by the injection momentum, pp cosOq. 



Since we are only interested in paths that have an 
intersection as small crossing angle, e, we can restrict 
the probabilities inside the integral to such paths by 
defining Ri ( R 2 ) as the phase-space position for the 
first (second) visit to the crossing occurring at time t\ 
(i 2 ). We can then write dR 2 = uf.sinediidt 2 de and set 
R 2 = (r 2 , fa) = (i"i, <f>\ ± e). Now we note that the loop 
cannot close before the path segments have diverged to 
a distance of order the system size L. When closer than 
this the two paths leaving an encounter have hyperbolic 
relative dynamics, and the probability of forming a loop 
is zero. Hence the duration of the loop must exceed 

T i (e)=A- 1 ln[e- 2 ], (26) 

This means that the probability that a path starting at 
Yo crosses itself at an angle ±e and then goes to lead to, 
multiplied by its injection momentum p-p cos0o> is 

roc pt—T\y/2 rtz—TL 

J m (Y ,e) = 2v F sine / dt dt 2 / dti 

JT L +T W JT l +T w /2 JT w /2 

x f dY [ dRi P(Y,R 2 ;<-i 2 ) 

Jm JC 

xP(R 2 ,R i; i 2 -h) P(Ri,Y ;ti), (27) 

where the Y-integral is over the cross-section of the mth 
lead, and Tw,Tl are shorthand for IV (e), Pi( e )- 

To get g wl ° we sum only contributions where 7I crosses 
itself, we then take twice the real part of this result to 
include the contributions where 7I avoids the crossing 
(and hence 72 crosses itself). Thus 

9tl = (^T 1 / dY /" 00 deRe[e i ^/ R ](7 m (Yo 1 e)>. 

J mo J 

(28) 

To average 7 m (Yo,e), we have to consider the average 
behaviour of the Ps. Within 7\y(e)/2 of the crossing the 
two legs of a self-intersecting path have a joint survival 
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I Usual weak-localization (wlO) 




Sucessful coherent-backscattering (cbsO) 




loop 



r 



Failed coherent-backscattering (cbsl) 




Failed coherent-backscattering (cbs2) 




FIG. 4: Weak-localization contributions to conductance (con- 
tributions to conductance of order iV ) when there are tunnel- 
barriers on the leads. All contributions except the second 
are for any mo and m, while the second (cbsO) is only for 
m = mo. The vertical cross-hatching indicate the regions in 
which the paths are correlated with each other, so their sur- 
vival probability is given by Eq. (|13|l not Eq. (|10|l . The en- 
counter occurs at the centre of the correlated region, as shown 
in Fig. [3] Details of the correlated regions for successful and 
failed coherent-backscattering are shown in Fig. [S] 



probability given by Eq. (|T3|) , since they are exploring 
the same region of phase-space. Elsewhere the paths' 
survival probability is given by Eq. (|10p . because there 
they are exploring different region of phase space. As a 
result the survival probability of a self intersecting path 
of length t > T W + T L is exp[-(£-27V)/<r m -T w /t D2 \. 
Since T m ' < 2r, 



Dl ! 



self-intersecting paths have an en- 
hanced survival probability compared to non-crossing 
paths of the same length. When the tunnel-barriers are 
absent, we have td 2 = Tdi and the survival probability is 
cxp[— (t — Tw)/tdi], as in Refs. 18ll20| and as first noted 
in Ref. Outside the correlated region, the legs can 



escape independently at any time through either lead. It 
is natural to assume that the probability density for the 
path to go to a given point in phase-space is uniform. In 
this case, the probability density for leg 1 (including the 
tunnelling into the cavity from lead mo) gives 

(P(Ri,Y ;ti)> (29) 
exp[-(ft - T w /2)/r m - 7V/(4t D2 )] 



p mo PF COS fn- 



the loop's probability density is 



2ttA 



(P(R 2 ,R i; t 2 -ft)) (30) 
exp[-[t 2 -ft_- T w /2]/t b1 - T w /(2t D2 )} 
2ttA 

and the probability density for leg 2 (including tunnelling 
into lead m) is 

<P(Y,R 2 ;i-ft)) (31) 
cos(9exp[-(t -ta - T w /2)/r ul - T w /(4t D2 )] 



>W„ 



Note that all the above probabilities are conditional on 
the fact the path has an encounter, so that the pair of 
path segments within Tyy /2 of that encounter have a joint 
survival probability given by Eq. (|13p (for convenience we 
divide that joint survival probability equally between the 
path segments). Thus we find that 

<P(Y, R 2 ; t - t 2 ) P(R 2 ,Ri; t 2 - h) P(R 1; Y ; ft)) 
PrriQ Prn Pf cos 9 cos 9q 
Z (2ttA) 2 2r D iEm' Pm>W m , 

x exp[-(jt - 2T w )/t u1 - T w /t U2 }, (32) 

so that (/ m (Y ,e)) becomes 

Pm pm (vf tb i ) 2 N m p F sin e cos O 



</ m (Y ,e)) = 



x exp[-Tvi//T D2 - (T L - T w )/tdi}. 



(33) 

We insert this into Eq. (|28|) . The integral over e is dom- 
inated by contributions with e <C 1, so that we write 
sin e ~ e and push the upper bound for the e-integration 
to infinity. Then the integral over e isi^, 



Re 



ede (Le/W) 



2/(Ar D2 ) 



cxp 



Xh 

2Ei 



-Rc 



jl+(ATD2) _ 



Xh 

r[i+(AT D2 r 1 



XhL 2 \ (Xrv2rl 



E F W 2 J 



2mw|,T D2 



cxphr^/TD^ + OpTD,)- 1 ], 



(34) 



where to get the second line we expanded to leading order 
in (Atd 2 ) _1 , in this situation the Euler Gamma-function, 
r[l + (Atd 2 )~ 1 ] — 1. Note that when we neglect all 
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C'[(Atd2)~ 1 ] corrections this must include all order-one 
terms in the logarithm of the Ehrenfest time, since they 
can only lead to 0[(Atd2) _1 ] corrections to the above 
expression. Thus since A ~ vp/L, we are justified in 
writing t° p = A" 1 ln[E F W 2 /(\HL 2 )]. 

Note that the second term in the exponent of Eq. (|33|) 
is independent of e, thus only td2 enters the e-integral. 
The Yo-integral generates a factor of 2W mo . We can 
write N m = (T7h)~ 1 p F W m and {J2 m > Pm'^m') -1 = 
(m>l) -1 ?iTDi. Thus the usual weak localization contri- 
bution is 

W 10 PmPm N m N mo Tpi 

9mm ° ' [E m ' Pm'N m ,}* T D2 

x exp[-r° p /r D2 - (r£ - r° p )/r D1 ], (35) 

where tdi,td2 are given in Eqs. (|llll4p . When the 
tunnel-barriers are transparent (p m = 1 for all m) , this is 
the only contribution to weak- localization for m ^ mo; in 
this case TD2 = Tdi and we get the expected resul t 16 i 18 ' 20 . 

Note that in the universal limit (Yg P — > 0) this contri- 
bution to conductance behaves like the naive argument 
given in Section Q] would predict. It goes like the weak- 
localization contribution to conductance for a cavity with 
no tunnel-barriers and each lead having N[ modes where 
N[ = PiN iy (with TDi/im being just a numerical prcfac- 
tor with 1 < tdi/td 2 < 2). Thus in the limit of opaque 
barriers (p m — > for all m), this contribution does not 
vanish for fixed Drude conductance (p m N rn remains con- 
stant for all m). Thus the origin of the suppression of 
weak-localization in the opaque limit is elsewhere. 



B. Successful coherent-backscattering. 

The coherent-backscattering contribution to trans- 
port has long been considered in trajectory-based 
semiclassics^i^. However until recently the approxima- 
tions were a little simplistic, and failed to capture the 
full nature of this contribution (in particular its Ehren- 
fest time dependence) . This was addressed in Ref. [HGJ] 
in the absence of tunnel-barriers. Here we perform a cal- 
culation similar to Ref. fl8| , while adding tunnel-barriers 
on the leads. The coherent-backscattering contributions 
in the presence of these barriers are shown in Fig. [4] 
with Fig. [5] showing the paths in the correlated region in 
more detail. We treat these contributions separately from 
those in the previous section because injection and exit 
positions and momenta are correlated. These contribu- 
tions are exactly those that were ignored in the previous 
section, paths where one/both legs escape within 7V/2 
of the encounter. However it is convenient to parameter- 
ize these contributions in terms of (roj_,po±) instead of 
(Mi)- 

Here we consider the successful coherent- 
backscattering (cbsO), the second contribution shown 
in Fig. [5] This successful contribution gives the usual 
coherent-backscattering when the tunnel-barriers are 



Sucessful coherent-backscattering (cbsO) 




Failed coherent-backscattering (cbsl and cbs2) 




FIG. 5: Details of the paths that contribute to successful and 
failed coherent-backscattering (cbsO, cbsl and cbs2). Path 
7I (solid black line) start on the cross-section of lead mo at 
position ?/o with momentum angle Go and returns at y with 
momentum angle 9. For the successful contribution it then 
tunnels into the lead, for the failed contribution it is reflected 
off the barrier and remains in the chaotic system. All paths 
are drawn in the basis parallel and perpendicular to 7I at 
injection, the initial position and momentum of path 7I at 
exit are r 0± = (y - y) cos 6 , r || = (yo - y) sin 9 and p 0± ~ 
-Pf(0-0o). 



absent (p m = 1 for all m) . From Ref. [HI , we have that 

S26 - Su = pft \\ + ^m\rl ± , (36) 
S 2 a-Si b = -pFr \\ +Poj_n>_L + \m\rl ± , (37) 

where (?"o_l,Po-l) are the position and momentum of path 
segment 7I6 relative to path segment 7la, shown in 
Fig. Thus the total action difference between these 
two paths is 

<5Scbso = {po± + m\r 0± )r j_ . (38) 

The successful coherent-backscattering contribution to 
the reflection is 

<4m° = S mmo (2nh)- 2 [ dYodY [°°dt 

J mo J 

x(P(Y,Y ;t)) Rc[e i5S< = bs0 ^], (39) 

where 5 mmo is a Kroneckcr (5-function. 

To perform the average we define Tyy(ra±,po±) an< A 
T^(rox,Po±) as the time between touching the tunnel 
barrier and the perpendicular distance between 71a and 
7I6 becoming W and L, respectively. For times less 
than Tyy{ro.LiPo±)i the path segments (71a and 7I6) are 
paired thus their joint survival probability is given by 
Eq. (fl3|) . For times longer than this the path segments 
escape independently, so each has a survival probability 
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given by Eq. (fTO]) . For g^mo we consider only those paths 
that form a closed loop, however they cannot close until 
the two path segments are of order L apart. This means 
that the i-integral in eq. (j3"9"|) must have a lower cut-off 
at 2T^(roxiPox)i hence we have 



dY / dt(P(Y,Y ;i)) 



2t; 



Pm PF COS 9 



x exp[-T^/r D2 



Pm'N„ 

- 2(n - 



where T' L w are shorthand for T' T 



{po± + m\ro±_) we estimat* 



,33 



L,W 



Tw)/t d1 ], 
iXW 



|po_L + m,Ar j 

mXL 
\poi_ + mXr j 



(40) 
For small 

. (41) 
. (42) 



Tw(r _L,p j_) - A 1 In 
T' L {ro±,Po±) ^ A -1 In 

Thus (T' L - T{y) 
deed 2(T[ - T> w ) 
above expressions into g^^ , write^ pFCOS# dYo = 
dyod(pF sin 0o ) = droj_dpo_L, and then make the substi- 
tution po = po± + m\ra±. We evaluate the ro^— integral 
over a range of order W mo . Then, taking the limits on 
the resulting po-integral to infinity, it takes the form 



is independent of (roj_,po±), in- 
= (rg 1 — Tg P ). We substitute the 



2hsm(p W/h) 
dp - 



Po 



l/(Ar D2 



4fi 



(mAPF) 1 /(Ar D2 ) 



Im 



mXW 

dp a p Q 1+1/(XTD2) eMmW/h] 



o 



contributions because they do not exist in the absence 
of the tunnel-barriers. They are part of the coherent - 
backscattering peak for the chaotic system, however they 
arc reflected back into the cavity by the tunnel-barriers. 
They can then contribute to either transmission or re- 
flection. 

The action difference for coherent-backscattering paths 
that reflect at the barrier is slightly different from that 
for the paths that transmit. This is because the path 
segments 7I6 and 726 Converge at infinity rather than 
at the lead (see Fig. [5]). We split path segment jib into 
7I6' before the failure to tunnel and 7I6" after the failure 
to tunnel. We do the same for path segment 726. Then 
the action difference is the sum of three terms. The first 
two come from before the failed tunnelling; they arc 

S2V - Si a = p F r || + |mA[r ±/2 - p Q j_/(2mX)} 2 .(45) 

and (S*2a — Sib') given by eq. (|3"T)) . The third term is 
the action difference of path segments after the failure to 
tunnel, 

Szb" -Sun = -Po±[ro±+Po±/{mX)] 

+±mX[r 0± +p ±/(mX)} 2 . (46) 

Hence the total action difference between the two paths 

SS chsl = SS chsQ + Po7(4mA) (47) 

where we write po = pox + mXr^x and hence SS c b s o = 
r o±Po, as in Section IV Bl As before we evaluate the inte- 
gral over roj_ first, getting an integral for po of the form 



2irh 



1XW 2 



l/(Ar D2 ) 



(43) 



To evaluate the integral we wrote it in terms of an 
Eulcr T-function. Note that [ft/ '(mXW 2 )} 1 ^^ ^ ~ 
exp[— Tg P /rD2] once we have dropped (D(l)-terms inside 
logarithm. The result is that 



cbsO 
ymmo 



Pmg Nmg 



x exp 



-r° p /rD 2 -(r^-r° p )/r D1 ]. (44) 



Note that the successful coherent-backscattering contri- 
bution has exactly the same exponential dependence on 
Tg P , Tg as the successful weak-localization. This must be 
the case if the theory is to preserve unitarity (conserve 
current), however since wlO and cbsO were calculated sep- 
arately using different methods, this acts as a check on 
our algebra. 



C. Failed coherent-backscattering. 

Here wc consider the failed coherent-backscattering 
(cbsl and cbs2), see Fig. @J:,d. We call them failed 



2ftsin(p W/ft) 
dp ; exp 



27rft 



h 



4mAft 



Po 



iXW 



1/(Atd 2 ) 



1XW 2 



l/(Ar D2 ) 



(48) 



To get this result we noted that po < h/W in the re- 
gion which dominates the integral, so we approximated 
cxp[ip§(4mAft)~ 1 ] = 1, after this the integral is identical 
to Eq. (PI). 

For the first failed backscattering contribution (cbsl) 
the bulk of the derivation is identical to that of the suc- 
cessful coherent-backscattering (cbsO) above. The dif- 
ference here is that at the point where the path re- 
turns to the lead it is reflected from the barrier (with 
probability 1 — p mo ) instead of being transmitted (with 
probability p mo ). The path then remains in the cav- 
ity, and will eventually escape through an arbitrary lead, 
with the probability of escape through lead m being 
E m ' Pm< N m i] p m N m . Thus we can conclude that 



cbsl 



cbsO x 



1 



p m N„ 



Pmo Xw Pm'N„ 
PmPrno 

(1 - p mo )N m N mo 



E, 



iN n 



x cxp[-T° p /r D2 - (t£ t° p )/t d1 ] (49) 
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The second failed cohcrcnt-backscattering contribution 
(cbs2) can be immediately written down by swapping all 
subscripts mo with to, hence 



cbs2 
"mmo 



x cxp[-r°P/r D2 - (rg - r°P)/r D1 ] (50) 



For to = too, one must confirm that cbsl and cbs2 
arc not double-counting the same contribution. One can 
see they form different contributions from the fact that 
one starts with a leg (two paths identical) and ends with 
a loop (two paths non-identical), while the other starts 
with the loop and ends with a leg. 



D. The total weak-localization correction with 
tunnel-barriers. 

We sum the contributions calculated in the the preced- 
ing subsections, Eqs. (|35I44I49I50|) . and get 



wl 



E, 



>N„ 



x exp[-r° p /r D2 - (r£ - r° p )/r D1 ], (51) 



where we define 

Anmg = Pm ~t~ < 



Em' Pm'Nr, 



m' 

(52) 



with S mmo being a Kronecker <5-function. Note that A 
goes like (p x A^ ). Thus in the limit p — > with con- 
stant pN, we see that .4^0 and the weak-localization 
contribution to conductance goes to zero. 

As a check on our result we verify that Em 9mm = 0> 
this means that we have not violated unitarity (or current 
conservation). 

Finally we consider the special case of two leads (L and 
R) with barriers with tunnel probability pl.r and number 
of modes ATl,r, then we find that Eq. ([STj) reduces to 



p 2 L p 2 R N L N R (N L + N n ) 
(plN l + p R N R ) 3 
x exp[-r° p /r D2 - (r£ - r° p )/r D1 ] (53) 



with Sgll 



VI. SHOT-NOISE FOR A CHAOTIC SYSTEM 
WITH TUNNEL-BARRIERS. 

Now we turn to the zero-frequency shot noise power, S, 
for quantum chaotic systems. This intrinsically quantum 
part of the fluctuations of a non-equilibrium electronic 
current often contains information on the system that 
cannot be obtained through conductance measurements. 



We give our results in terms of the Fano factor F = S/S p , 
which is the ratio of S to the Poissonian noise S p = 2e(I) 
that would be generated by a current flow of uncorrelated 
particles. Alternatively, one can think of the Fano factor 
as the ratio of the noise to the average current, writ- 
ten in convenient dimcnsionlcss units. According to the 
scattering theory of transport one has^l 



p = Tr[S mmo §mm ] - Tr[St 



It 



'mm a mm ^ m m Q mm 



Tl [Smmo ^mmo J 



.(54) 



We use the trajectory-based method developed in 
Ref. to evaluate the contributions to this quantity 
to lowest order in A^" 1 , where N is the number of lead 
modes. We thereby neglect all the weak-localization-type 
corrections to the noise calculated in Ref. [|[ . At this or- 
der, all contributions are listed in Figs. [6] and [Jj The 
contributions fall naturally into two classes. The first 
class, shown in Fig. [6l involve classical paths (path 1 
and 3) which are correlated, and one or both paths es- 
cape before their flow under the cavity dynamics makes 
them become uncorrelated. The second class, shown in 
Fig. [Jj involve correlatcd-paths which become uncorre- 
lated under the cavity's classical dynamics before any 
paths escape. As before, correlated means that the paths 
are almost parallel and within W of each other. 

The denominator and the first term in the nu- 
merator of Eq. ([Ml) arc equal to the Drude conduc- 
tance, and are hence given by Eq. (|22[) with to ^ 
too- Thus to find the Fano factor we must evaluate 
Tr[Sj„ mo S mlI1(1 §!„ mo S mmo ]. As with weak-localization wc 
approximate ^ n {y'\n){n\y) « 5{y' — y), then it becomes 
a sum over four paths, 7I from j/01 to yi, 72 from j/03 
to yi, 73 from yo3 to 2/3 and 74 from j/01 to j/3, where 
yoijyo3 are on lead toq and 2/1,2/3 are on lead to. Hence 



T r [§>mm ^mmo^mmo ^ 
1 



(2^)"JL 



d2/oidj/o3 / d2/id2/ 3 

R 



(55) 



x A* l4 A l3 A* 2 A 7l exp[iSS/h] 

7l,---74 

where A 7 is given by Eq. (fT?)) . and SS = S^i~ S , 7 2 + 5 , 7 3 — 
yS 7 4 (we have absorbed all Maslov indices into the actions 
S-yi). The dominant contributions that survive averaging 
over energy or cavity shape are those for which the fluctu- 
ations of SS/h are minimal. They are shown in Figs. [6l7l 
Their paths are in pairs almost everywhere except in 
the vicinity of encounters. Going through an encounter, 
two of the four paths cross each other, while the other 
two avoid the crossing. They remain in pairs, though 
the pairing switches, e.g. from (71574) and (72573) to 
(71572) and (73574). Paths are always close enough to 
their partner that their stability is the same. Thus, for 
all pairings 



E 

71,. . .74 



71,73 



1-4 



-1 



(56) 
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(8)2*0 



path 1 




path 3 





FIG. 6: The set of contributions to tr[Sj rlmo S mmo Sj! nmo S mmo ] which we call D la _, ■ ■ ■ , D ri . These are the contributions which 
do not vanish for infinite Ehrenfest time. Here we show only the tunnel-barriers on leads mo and m as shaded rectangles, a 
path which crosses the barrier on lead m has succeeded in tunnelling out of the cavity into the lead. The contributions are 
made up of four classical paths, here we show only two of the paths (1 and 3). The other two paths (2 and 4) look the same as 
the paths shown, except that they cross at the centre of the correlated region (indicated by the vertical crosshatching). In the 
correlated region the paths have the same topology as shown in Fig. [3] if we replace the labels in the manner given in Eq. (|58|l . 
Thus path 4 is paired with path 1 at lead mo but paired with path 3 at lead m (and vice- versa for path 2). The noise in these 
contribution is purely due to the stochastic nature of scattering at the tunnel-barriers, if they were absent these contributions 
would be noiseless. The calculation will show that all these contributions go like (1 — e~ T B / TD2 ). 



Then the sum over all paths 7 from yo to y is 
given by Eq. (|20[) As with weak-localization, we define 
P(Y, Yo; t)SyS95t as the product of the momentum along 
the injection lead, pf cos^o, and the classical probability 
to go from an initial position and angle Y = (yo, ) to 
within (Sy, 56) of Y in a time within St of t. We remind 
the reader that this probability is to go from the injec- 
tion lead to the exit lead, and this includes the tunnelling 
probability at each barrier. We now use Eqs. (|55|) . ([55]) 
and (j2"0)) to analyze the contributions in Figs. 16171 All 
contributions can be written as 



make the replacements 



A = 



1 



dY 01 dY 03 / dYj dY 3 dh dt 



(27Tft) 2 jL ,R 

x (P(Y 1 ,Y 01 ;t 1 )P(Y3,Y 03 ;t 3 )) cxp[iS S D t / h] , (57) 

where the subscripts 1,3 indicate paths 1 and 3 respec- 
tively. When evaluating Eq. ([57]) the joint exit probabil- 
ity for two crossing paths has to be computed. 



Noise contributions which do not vanish for 
infinite Ehrenfest time. 



To evaluate all the contributions in Fig. [HI we note that 
the paths never become uncorrelated under the classical 
dynamics; they only escape in an uncorrelated manner if 
one path tunnels while the other is reflected. In this case 
the details of the encounter are as given in Fig. [5J if we 



path segment 71a — > path 7I, 
path segment 72a — ► path 72, 
path segment 7I& — > path 73, 
path segment j2b — > path 74. 



(58) 



Thus the action difference between the paths can be eval- 
uated in a manner equivalent to cohcrcnt-backscattcring. 
For contributions where both paths escape as a pair, the 
action difference is given by SSoi = ^cbso where 8S c \> s o 
is given in Eq. (|38p . For contributions where the pair is 
broken by one of the paths escaping, the action difference 
is given by SSoi = ^cbsi with SS c h s i given in Eq. (|47[) . 
Here, as in Section IV C[ the difference between 5S c t, s i 
and <5Scbs0 nas 110 effect on the final result, so one could 
just use SS c bsO for all contributions. 

For the contribution in Fig. [5^, the paths paired at lead 
mo remain paired at lead m, thus the length of the paired 
paths must be less than T^ v (ro± 1 po±) given in Eq. (|41[) . 
Thus we see that 



dY!dY 3 



dhdts (P(Yj , Y01 ; ii)P(Y 3 , Y 03 ; t 3 )> 



■(l-exp[-T^/r D2 ]). (59) 



p 2 m N m pl cos 6>qi cos 6>q 

Em' P™'( 2 - Pm') N m 



where is the function of (ro_L,poJ_) given in Eq. (|4ip . 
Note that the denominator comes from the fact we 
are considering the survival probability for the corre- 
lated paths; thus the probability that the correlation 
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is destroyed by paths escaping into a lead during the 
time t to t + 5t is P^it) x St/r^, where P2{t) is given 
by Eq. [fl3|). We inse rt Eq. ((59j) into Eq. [[ST]), then 
just as in Section IV Bl we change integration variables 
using pFCOS#o3dYo3 = dro±dpo±, arid define po = 
poi_ + m\ro±. In the regime of interest T{ v (ro±,po±) — 
A -1 ln[mAHY|po|]- Evaluating the integral over roj_ 
leaves a po _m tegral which we cast as an Euler T-function, 
just as for cohcrent-backscattering in Section IV Bl To 
lowest order in (At^) -1 we find the contribution to 



Lmo^mmo^mmn^mrnJ shown ill Fig. [(Jl is 



P 2 ma P 2 m N m N m {l - exp[-rg p /r D2 ]) 

Em' Pm'( 2 - Pm>)N m > 



(60) 



Now we note that each contribution in Fig. [6] is of 
a similar form to Z?i a . For example, -Dib and D\ c are 
like -Di a with the exception that a path is reflected off 
lead m and then returns to lead m. The result of the 
integral over (tq±,pqj_) is basically unchanged when we 
replace the action difference in Eq. with Eq. (|37]l. 
see Section IVCl Thus we get each of these contributions 
by simply multiplying D ia by 



1 - P,i 

Pm 



p m N„ 



Em' Pm'N„ 



(1 - Pm )N„ 
Em' Pm'N m 



(61) 



where after reflection the path that remains in the cavity 
evolves alone, so its survival is governed by Eq. (TIT))) . We 
can make the same argument for paths which enter the 
cavity from lead nig at different times, but in such a 
way that the second enters the cavity at a moment when 
by chance the first is reflecting off barrier mo, in such a 
way that the paths form a pair. To make the argument 
we need simply reverse the direction of the paths, and 
we return to the situation discussed above Eq. (jSTj) with 
m replaced by rriQ. Thus we find that the sum of all 
contributions in Fig. [6] is 



D, = 



P 2 m p 2 m N mo N m (l 



exp[ 



Em' Pm'(2 ~ p m >)N„ 

| 2{l- Pmo )N ma 

Em' Pm'N m < 



(62) 



2(l-p m )iV„ 

Em' Pm'N m - 



B. Noise contributions which vanish for infinite 
Ehrenfest time. 

To evaluate all the contributions in Fig. [7] except 
Fig. [7g, we note that the path pairs are always corre- 
lated at one lead; at that lead they only escape in an 
uncorrelated manner if one path tunnels while the other 
is reflected. In this case the details of the encounter are 
as given in Fig. [5l if we make the replacements given 
in Eq. (|58p . These contributions are extremely similar 
to those discussed in Section [VI Al the only difference is 
that the paths become more than the lead width apart at 



at least one lead. This occurs because the paired paths 
survive for a time longer than T^,. Thus for the contri- 
bution in Fig. [7^, we have 

I dYidY 3 / dt 1 dt 3 (P(Y 1 ,Y 01 ;t 1 )P(Y 3 ,Yo3^3)) 

Jrn JT^ 

_ Pm N mPF cos 001 COS 9 03 



(Em' Pm'N m ,) 2 



exp[-T^/r D 2] 



(63) 



Since the paths separate inside the cavity under the clas- 
sical dynamics, their escape is given by the single path 
escape. However the factor of exp[— T^/t^] is due to 
the fact that the correlated paths must first survive un- 
til time T^y(ro±,Po±), and during this time the survival 
probability is given by Eq. (fT3"|) . Compare this with the 
equivalent contribution, Eq. (|59|) . for pairs which escape 
on a time less than T^y. The integral over (roj_,poJ_) can 
be evaluated in exactly the same manner as D\ a (and in 
the same manner as the coherent-backscattering in Sec- 
tion |VB|. Then we find that 

n p 2 m p 2 m N mo Nlexp[-T° P /T D2 } 
~ (Em'P-^') 2 ' 

We then note that contributions D 2 b an d D 2c equal D 2& 
multiplied by (1 - Pm )/(Em/ Pm'N m >), thus 



£>2a + D 2h + D 2c 



P 2 mo P 2 m N m iV m exp[-T° p /r D2 ] 



X 1 + 



(Em' p m >N m ,y 

2(1- Pmo )N ma 
Em' Pm'N m i 



(65) 



The contributions D 2( \, D 2e and D 2 i (shown in 
Fig. [7ji,e,f) are the same as contributions D 2a , D 2 \, and 
D 2c (shown in Fig. [7k-,b,c) with mo <-> m. Thus 



D 



2d 



D 



2 c 



P 2 mo P 2 mN^ N n 



exp[ 



op 



(Em' Pm>N m ,y 

2(1 - Pm )N, 

Em' Pm'N„ 



(66) 



This leaves us to evaluate D 2g , the contribution shown 
in Fig. [7^, this is different from all other noise contri- 
butions because the paths are not correlated at escape. 
Thus we evaluate this contribution in a manner simi- 
lar to the weak- localization contribution in Section [V Al 
We use the method developed by Richter and Siebcr 5 , 
while taking into account that paths in the same re- 
gion of phase-space have escape probabilities given by 
Eq. ([TBI. This method was first applied to shot-noise 
in Ref. [atl7| in the absence of tunnel-barriers. Here we 
follow Ref. [13, and we write the action difference as in 
Eq. (|2~4"]) . where the crossing angle, e, is shown in Fig. [7g. 
We write 

P(Yh Yoi; U) = / dRjP(Yi, R,;; t\ — t^)P(Rj, Yoi; tj) , 
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path 3 
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FIG. 7: The set of contributions to tr[S} rlmo S mmo Sj! nmo S mmo ] which we call D 2a , • ■ ■ , D 2s . These are the contributions which 
vanish for infinite Ehrenfest time. The contributions are drawn in the same manner as Fig. [6] with only paths 1 and 3 shown. 
As before paths 2 and 4 are equivalent to 1 and 3 except that they cross at the centre of the correlated region (indicated by 
vertical crosshatching) . For finite Ehrenfest times, the paths must escape the correlated region (vertical crosshatching) before 
either tunnels out of the cavity (on at least one side of the encounter). The calculation will show that all these contributions 

op / 

go like e _T E /™ . 



where P is the probability for the classical path to exist 
(not multiplied by the injection momentum), and R^ is a 
point in the system's phase-space fa, fa) visited at time 
t' i: with fa giving the direction of the momentum. Wc 
choose Ri and R3 as the points at which the paths cross, 
so R 3 = (n, fa ± e) and dR 3 = Vp sined^d^de. 

To get D 2 g we sum only contributions where 7I crosses 
73, we then take twice the real part of this result to 
include the contributions where 7I and 73 avoid crossing 
(and hence 72 and 74 cross). Thus 

D 2g = 2(27rfr)- 2 ^dYoidY 03 

x I" deRe[e iSS °i/ H ](l(Y 01 ,Y 03 ;e)). (67) 
Jo 

where <5<S(di i s the same as SS^i in eq. (|24[) . The function 
/(Y01 , Y03; e) is related to the probability that 73 crosses 
7I at angle ±e. Its average is independent of Y i,o3, so 
(T(Y i, Y 03 ; e)} = (/(e)). Injections/escapes are more 
than T\y(e)/2 from the crossing, so 



(/(e)) = 2v\ sine J dY^Y z J dRi 



ptt-T/2 poo pt 3 -T/2 

I dh / dt[ / dt 3 / dt' 3 

IT JT/2 JT JT/2 



(P(Yi,Ri;ti-ti)P(Ri,Y i;ti) 
P(Y 3 ,R 3 ;< 3 - < 3 )P(R 3 , Y 03 ;t 3 )) , (68) 



where TV is the function of e given in Eq. (|2"3"|) . We 
next note that within TV/2 of the crossing, paths 7I 



and 73 are so close to each other that their joint survival 
probability is given by Eq. (|13p . Elsewhere 71, 73 escape 
independently through either lead at any time, hence 



, _ p F r D1 P mo Pm N m cos 6»oi cos (9 03 sine 



xfvm 



(Em' Pm'N m ,f 



-T w (e)/T 02 



(69) 



where we used N m ' = {'nh)~ 1 pFW m ' , and assumed that 
the probability that 73 is at R 3 at time t' 3 in a system 
of area A is (2ttA)~ 1 = m[2irhTpj Em' Pm'-^m'] -1 - Then 
the Y i i03 -integral in Eq. J67]) gives (2W mo ) 2 - The in- 
tegral over e is the same as given in eq. (f3~4"|) for weak- 
localization. Thus 



Do 



Pm PmN mo Nm Em' Pm> P ~ P m ')N n 
(E m ,Pm'N m ,)* 

x exp[-r° p /T D2 ] ■ 



(70) 



Summing all the contributions in Fig [SJ given in 
Eqs. (|65l66l70p . we get 



Do = 



Pm Pm N m N m exp[-Tg p /-r D2 ] 
(Em' Pm'N m ,f 

2(2 - pm - pm)Nm N„ 



N mo + N m + 



E/m' Pm'N m r 
N mo N m J^m' Pm' (2 ~ Pm')N„ 
(Em' Pm'Nm') 



(71) 
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C. Total shot-noise for arbitrary Ehrenfest time. 



Thus we find that the Fano factor is given by 



1 - 



Prao Pm mo n 



(72) 



where D\ and D% are given by Eqs. (|63I71I) . 

Let us first consider this result for transparent barriers 
(Pm' = 1 for all m'), in this case it simplifies to 



F = 



1 



N m N mo N m 



(En.'-Nm') S 



ex P[- T E P / T D] 



(73) 



where r D 1 = (toL)" 1 £^ =1 W m . As noted in Ref. g||l7| 
for two leads without tunnel-barriers, the contributions 
to Tr[§J nmo S mmo Sj )imo S mmo ] which enter or exit the cav- 
ity in a correlated manner cancel the the contribution of 
Tr[§J nmo § m „ lo ]. In that case the noise is given by the 
contribution which enters and exits the cavity in an un- 
corrected manner (Fig. [6g). Here we see this is only 
the case for two leads; the first term in Eq. ([73|) comes 
from Tr[SJ„ m S mmo ], while the second term comes from 
those contributions to Tr[8 mTOo S mmo S}„ mo S mmo ] which 
enter or exit the cavity in a correlated manner; in general 
these two do not cancel each other. However we see that 
for an arbitrary number of leads without tunnel-barriers 
the paths shorter than the Ehrenfest time are noiseless 
(just as for two leads^). 

Now let us consider eq. (f72| in the limit of opaque 
barriers (p m i — > for all to'), wc see from Eq. (|72p that 
the Fano factor is 



F = 1 



2Pm a PmN ma N„ 



(74) 



which is completely independent of the Ehrenfest time. 
In fact, when we expand Eq. (f72|) to 0[p], we see that 
the 0[/9]-term is also independent of the Ehrenfest time, 
thus 



F = 1 - 



N„ 



(E m < Pm'N m ,f 
Pm PmN mQ N m ][] m , p 2 m ,N m , 

{Y,m'P m >N m >Y 
PrrioPm 

(N ma + N m ) „ 



Em' Pm'N„ 



o[ P 2 



(75) 



where the 0[p 2 ]-term is dependent on the Ehrenfest time 
In the case of a symmetric two lead cavity, JVi = N 2 = 
N and pi = p 2 = P, the rather ugly result in Eq. ([72| 
greatly simplifies to 



F=^—^(l-c-< P /^ 



2 ~P c -rgP/r D2 



(76) 



Both terms in this expression increase monotonically as 
we reduce the transparency of the tunnel-barriers from 



p = 1 to p = 0, thus the presence of the tunnel-barriers 
always increases the quantum noise (for given Ehrenfest 
time). In the limit of transparent barriers (p = 1), the 
Fano- factor crosses over from the RMT value of 1/4 to 
zero as we increase the Ehrenfest time. While in the 
limit of opaque barriers (/) « 1) the Fano factor is 1/2, 
independent of the Ehrenfest time. 



D. Consistency check: alternative shot noise 
formula. 



Finally we note that the main difficulty in the above 
calculation is the identification of all contributions. Thus 
it is important to have an independent check that no 
contributions have been missed. In appendix [Cl we use 
the fact that the Fano factor can be written in terms of a 
product of transmission and reflection contributions, see 
Eq. (|C2[) . and rederive the result in Eq. (|72|) starting from 
that formula. The transmission/reflection contributions 
calculated there (shown in Figs. [8] and [9]) combine in a 
non-trivial manner to give the result in Eq. (|72jl . This 
can be seen by the fact that there is no obvious too <-> rn 
symmetry in the contributions in Appendix [C] (compare 
that with the contributions calculated above, where every 
contribution has a partner with to <-> to). Thus if wc 
miss a contribution in Figs. IBTTl and miss the equivalent 
contributions in Figs. 18191 it would be extremely unlikely 
that the two sets of contributions would sum to give the 
same result. 



VII. CONCLUDING REMARK. 

We have summarized the central physics discussed in 
this article in Section [Til an( l recommend that as a sum- 
mary of the contents of this article. Here we would like 
to re-iterate one technical point which makes the semi- 
classical calculations tractable. It is the fact that when 
writing probabilities to escape into a given lead wc count 
only successful attempts to escape (not situations where 
the particle hit the tunnel-barrier on a lead and was re- 
flected back into the cavity). This is very different from 
the equivalent RMT calculation^! in which each collision 
with the tunnel-barrier on a lead was explicitly taken 
into account, regardless of whether the particle tunnels 
or not. We feel unqualified to speculate if our approach 
could be directly applied to simplifying the RMT calcu- 
lations. However we are certain that the hand-waving 
arguments made in Section III Al would be an excellent 
guide to the intuition when performing involved RMT 
calculations. 
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APPENDIX B: APPROXIMATING SUM OVER 
LEAD MODES BY DIRAC 5-FUNCTION. 



APPENDIX A: PHASE OF TRANSMISSION AND 
REFLECTION AT A TUNNEL-BARRIER. 



Here we calculate the complex transmission and reflec- 
tion amplitudes, t and r, for a rectangular tunnel-barrier, 
in the limit that the barrier is thin and high (width I and 
height U) . We follow the standard procedure explained in 
most quantum mechanics textbooks. We include the de- 
tails here simply because most textbooks give the trans- 
mission and reflection probabilities, but not the phases. 

We consider a one-dimensional system with the par- 
ticle comes from the left. To the left of the tunnel- 
barrier the solutions of Schrodinger's equation have mo- 
mentum ±p, while to the right of the barrier the only 
non-zero solution has momentum +p. Inside the bar- 
rier the solutions decay/grow like exp[±kx/h], where 
k = [2mU — p 2 ] 1 / 2 . Solving the simultaneous equations 
that we find by matching the wavefunction and its deriva- 
tive at the two edges of the tunnel-barrier, we see that 
the transmission amplitude is 



exp[— ipl/h] 



cosh{kl/h) + Uk/p-p/k) smh{kl/h) 



while the reflection amplitude is 



(k/p + p/k) s'mh(kl/h) 



2cosh(kl/h) + Uk/p-p/k) smh(kl/h) 



(Al) 



(A2) 



where Z is the width of the barrier. To take the thin 
high barrier limit, we take U — » oo and I — > such that 
kl/h remains finite. We then adjust kl/H to give us the 
tunnelling probability, p = \t\ 2 , that we wish; in this limit 
pl/h = 0. 

We will need the transmission and reflection phases in 
Appendix [Cj because for _Dib we find that t and r do 
not appear in the combination \t\ 2 = p or \r\ 2 = 1 — p. 
Instead we need to calculate the combinations (t*r) 2 and 
(r*t) 2 . From the above results for t and r we see that 



[( e W/»t)* r ] : 



[{e ipl/R t)r*} 2 



-\t\ 2 \r\ 2 = -p{l-p) (A3) 



and hence in the thin barrier limit (pi -C K), we get 
(t*rf = (r*t) 2 = -p(l - p). 

We note that throughout this article, we assume that 
p = \t\ 2 is unchanged as we take the classical limit 
Xp/L — > 0. This means that we keep the ratio Xp/l 
fixed (and do not scale Z like the classical scales W, L) 
as we take the limit Xf/L — ► 0. For simplicity we call 
this the classical limit. Now, one might argue that the 
strict classical limit would involve taking Ap/Z — > 0, so 
there would be no tunnelling. However this limit is un- 
interesting because without tunnelling there would be no 
conduction at all. 



The approximation made above Eq. (|19p was intro- 
duced in Refs. 17lfT§ |. however the explanation given 
there was too brief to be helpful. The issue not dis- 
cussed there was the fact that exp[i<S 7 /7i] in Eq. (fTo)) 
oscillates as fast with y as (j\y) does. Thus we must 
deal explicitly with both fast oscillating terms at once. 
To do this we write path 72 (which goes to y ) in 
terms of equivalent path that goes to y, then S^iy') = 
Sj2(y)+PF sin #'(?/ — y) + 0[mX(y' — y) 2 ]. Hence when we 
write tr[§J Tlrno S mmo ] we find that it contains (y 1 — j/)-terms 
which oscillate fast (oscillate on a scale h/p-p = Af), these 
terms are J2 n (y' \ n ) ( n \y) exp[ipF sm9'(y' - y)/h]. 

Evaluating the sum for an ideal lead with N lead modes 
of the form (y\n) = [2/W) 1 ' 2 sm(i:yn/W), one finds that 



]T)(y'|n) (n\y) 

n 

sin[(z' -z)(N- 



1/2)] 



2Wsm[(z' - z)/2] 



(Bl) 

sin[(z / + z)(iV+l/2)] 
2VFsin[(z' + z)/2] 



where z = ny/W. This function is peaked at y' = y 
with width ~ Ap and height ~ Ap 1 , its integral over 
y' is one. Hence ^2 n {y'\n){n\y) ~ S(y' — y) for func- 
tions which varies slowly on the lengthscale of a Fermi 
wavelength, where S(y') is a Dirac (^-function. How- 
ever the crucial point for this derivation is that one can 
also show that ^2 n {y'\n){n\y)exp[ipFSm9'(y' — y)/h] ~ 
8(y' — y). To see this we consider the integral, / = 
J dy'J2n(y'\ n )( n \y) ex P[iPFSm9'(y'-y)/h}. Defining z = 
ir(y' ~ y)/W, we find that 



(27T) 



dz 

z 



a i[iV(l+sin0') + l/2]2 



-i[JV(l-sin6»') + l/2]2 



where we use the fact that the integral is dominated by 
z ~ iV _1 <C 1. The integrand is finite at z = even 
though the individual terms in the integrand diverge, so 
we can push the contour of integration infinitesimally 
into the lower half of the complex plane. To evaluate the 
integral over each term in the integrand, we note that 
the first (second) term converges in the upper (lower) 
half-plane. The first term's contour is deformed into the 
upper-half plane, but then encircles the pole at z = 
giving the contribution 2n. The second term's contour 
encircles nothing when pushed into the lower-half plane, 
so contributes nothing. Thus 1 = 1, independent of the 
prefactors in the exponents. The integral is dominated 
by z ~ iV" 1 (i.e. y ~ Ap), hence 

^2(y'\n)(n\y) cxp[ip F sm6'(y' -y)/h] ^ S(y' -y), 

n 

for functions which vary slowly on the scale of the Fermi 
wavelength. 
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(a) 



(b) m 



path 1 





(d) m Q path 1 



m 



(f) 




m T^m 




(g) m Q P at h 1 



m &m 



FIG. 8: Here we list all the contributions to S m /^ m tr[S] M , mo § m ' mo S] Tlmo S mmo ] which remain finite in the limit of infinite 
Ehrenfest time. The contributions are drawn in the same manner as Fig. [6] with only paths 1 and 3 shown. The noise in these 
contribution is purely due to the stochastic nature of scattering at the tunnel-barriers, if they were absent these contributions 
would be noiseless. The calculation will show that (a) and (b) are independent of t^ p , while all other contributions go like 

(l-e- T E P /-D2). 



APPENDIX C: CALCULATING SHOT NOISE 
FROM TRANSMISSION /NON-TRANSMISSION 
CORRELATIONS. 



The unitarity of the scattering matrix means that 



E 



\m' mo^m' mo 



= i. 



(Cl) 



It follows that we can write the formula for the Fano- 
factor in Eq. ([54]) as 



F = 



Hm'=£m T r [^m'm ^ m ' mo ^mmo ^"im ] 



Tr[S. 



(C2) 



mmo ^mmo J 



Thus we see that the Fano factor is given by correlations 
between transmitting contributions (those going to lead 
m) and non-transmitting contributions (those going to 
any lead except m). 

We now use the trajectory-method to calculate this di- 
rectly. We find no preservation of unitarity at the level 
of the individual path's contributions to Eq. ([54"]) and 
individual path's contributions to Eq. (|C2|) . it is only 
when we have summed all contributions to Eq. (|54j) and 
Eq. (|C2[) . that we find them to be equal. We use this as 
a check on the fact we have included all relevant contri- 
butions, and as a check on our algebra. 

To calculate these contributions we use exactly the 
same method as in Section I VII For a given contribu- 
tion, one can take a contribution with the same topology 



in Fig. [5] or Fig. and note that one of the legs goes 
to m! rather than m, with ml being summed over all 
leads except lead m. Thus here we simply give the re- 
sults for each contribution without discussing details of 
the calculation. The exception to this is the contribu- 
tion in Fig. [SJd which has no analogous contribution in 
Section IVTl 



The first two contributions in Fig. [5] arc independent 
of the Ehrenfest time, they arc 



Pmp ( 1 ~ Pmp ) Pm N mo N m 

2p m( M - Pm a )p m Nl N„ 
(Em' P m 'Nm'f 



(C3) 
(C4) 



Note that the negative sign in is due to the tunnelling 
and reflection phases at the tunnel-barrier on the moth 
lead. Here, unlike for all other contributions, the complex 
tunnelling amplitude, t, and reflection amplitude, r, do 
not appear as products of |r| 2 = p and \t\ 2 = (1 — p). 
Instead they appear in the combination (t*r) 2 or (tr*) 2 , 
thus we need the phases of t, r. These are calculated 
in Appendix [Al where we show that for the high, thin 
barriers considered in this article the phases generate an 
overall negative sign, so (t*r) 2 = (tr*) 2 = —p(l — p). 

All other terms in Fig. [5] go like (1 — exp[— t^F /t^])- 
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(c)m 



path 1 m 



m ^tn 




m 



(d) »„ ^.Pathl 




m 



(f) m „^pathl 




tri 



which vanish in the limit of infinite Ehrenfest 



FIG. 9: Here we list all the contributions to ^ mVm tr[S^, mt) S m / mo § mmo a mmoj 
time. The contributions are drawn in the same manner as Fig. [6] with only paths 1 and 3 shown. The calculation will show 
that all these contributions go like e _T E / T D2 Thus for t£, p /td2 ~* 0, the noise is given by the sum of these contributions and 
Fig. [9ji,b. 



The sum of the next three contributions in Fig. [5] gives 
Pm PmNm N m (l - exp[-T° p /r D2 ]) 



The other three contributions in Fig. [5] arc 



-Dl(c+d+c) 



(Em' Pm'( 2 - Pm')N m ')(£, m , p m 'N m >) 

X y^/Q ro '(l - Pm>)N m i - p m (l - p m )N, 
\ m' 

2(1 -p mo )N mo 



D 2d = 



pLoPmN^Nrn exp [-t° p /t D2 ] 



X 1 + 



Em' Pm'N m ' 

The sum of the final three contributions in Fig. [5] gives 
P 2 mo P m N mo N m {\ - exp[-TE P /T D2 ]) 



(C5) 



(Em' Prn>N ml f 

^ Pm< (1 - P m ' )N rr , 
\ m' 



(C8) 



Pm(l - Pm)N v 



D 1 



(f+g+h) 



(Em' Pm'( 2 " Pm') N m')(Y, m ' P™' N m >) 
- Pm(l - Pm)N., 

2(l- Pmo )N mo 



(1 - pm)^Pm'N„ 



X 1 + 



Em' Pm'-Nm' 



(C6) 



The sum of all the contributions shown in Fig. [5J gives 
the noise in the infinite Ehrenfest time limit (when all 
contributions in Fig. [5] are zero), and is in agreement 
with the equivalent result, Di, in Section [VTl 

We now turn to the contributions in Fig. [9l The sum 
of the first three contributions in Fig. [9] gives 



D 



2(a+b+c) 



(C7) 



P 2 mo p m N mo N m cxp[-rE P /T D2 ] 



x 1 



Em' Pm'N„ 

PmNm 
Em' Pm'N m 



2(1 -p mo )N mo 

Em' Pm'N m i 



Do 



pLoPmN^ Nrn exp[-T° p /T D2 ] 



(Em' Pm.N m ,f 

(1 - p m ) ^ p m >N m 
\ m' 



(C9) 



Pm(l - Pm)Nr, 



D 



21' 



P 2 mo PmN^ N m exp[-T° p /r D2 ] 



(Em' Pm>N m ,y 



(CIO) 



x 1 



2_u m i Pm'^m' J . 



The sum of all the contributions shown in Fig. [5] and the 
contributions shown in Fig. [8jt,b gives the noise in the 
zero Ehrenfest time limit (when all other contributions in 
Fig. [8j are zero) , and is in agreement with the equivalent 
result, £> 2 , in Section PVTl 



Substituting the sum of all contributions in Figs. [Hand 
[9] into the numerator of Eq. (|C2[) , we again arrive at the 
result given in Eq. (|72|) . This shows that our method 
respects the unitarity of the scattering matrix (thereby 
conserving current), and acts as a check on the algebra. 
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